Crossover from weak localization to weak antilocalization in a disordered microbridge 
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We calculate the weak localization correction in the double crossover to broken time-reversal and 
spin-rotational symmetry for a disordered microbridge or a short disordered wire using a scattering- 
matrix approach. Whereas the correction has universal limiting values in the three basic symmetry 
classes, the functional form of the magnetoconductance is affected by eventual non-homogeneities 
in the microbridge. 
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Interference of time-reversed paths causes a small neg- 
ative quantum correction to the conductance of a disor- 
dered metal termed the weak localization. L2uLa This cor- 
rection is suppressed by a time-reversal symmetry break- 
ing magnetic field, whereas in the presence of strong spin- 
orbit scattering, the sign of the correction is reversed™ 
In that case, the interference correction is known as weak 
antilocalization. 

In a wire geometry at zero temperature, the weak lo- 
calization correction takes a particularly simple and uni- 
versal form™ 



SG = 



2e 2 (/3-2) 
3/3/i ' 



(1) 



where the symmetry parameter j3 denotes the appro- 
priate symmetry class: In the presence of an applied 
magnetic field, (3 — 2, and without a magnetic field, 
(3 = 4 or 1 with or without strong spin-orbit scat- 
tering, respectively. Equation was obtained using 
random-matrix theory^*^ and diagrammatic perturba- 
tion theory*^ and is valid if the length L of the wire is 
much smaller than the localization length £ and the de- 
phasing length L$, but much larger than the mean free 
path I. The validity of Eq. (JIJ extends to the case when 
sample parameters are non- homogeneous, e.g., for wires 
of varying cross section, mean free path, or electron den- 
sity. 

For wires with weak spin-orbit scattering, a crossover 
between weak localization and weak antilocalization 
takes place when the spin-orbit scattering length l so be- 
comes comparable to L or (whichever is smaller) . Ex- 
perimentally, this crossover regime has been well studied 
in wires with length L 3> L^MtlZ^ In this regime, weak 
(anti) localization takes the form of a small correction to 
the conductivity of the wire, rather than of a correction 
to the conductance. Theoretically, the weak localization 
to weak antilocalization crossover in the regime L 3> L<f, 
has been considered in Refs. Il llrlnl using diagram- 
matic perturbation theory. The opposite regime L -C 
where the universal correction Q to the conductance G 
can be observed, would be relevant for relatively short 
high-purity metal wires^i or disordered microbridges. 



The goal of this paper is threefold: (i) to general- 
ize the random-matrix methods for quantum wires to 
the crossover between weak localization and weak an- 
tilocalization, thus extending the equivalence of the two 
methods to the interpolation between the three sym- 
metry classes, (ii) to find an explicit expression for 8G 
for L -c L^, and (iii) to extend the theory for the 
crossover regime to the case of non- homogeneous wires, 
for which the electron density, impurity concentration, 
or cross section varies along the sample. In this case, 
both the crossover scale and the functional form of 6G 
in the crossover are affected by non-homogeneities. The 
fact that the crossover scale, characterized by the spin- 
orbit length l so and the magnetic length Ih, is non- 
universal is well known, both for homogeneous and for 
non-homogeneous microbridges^ Our finding that the 
functional form of the crossover is affected by the non- 
homogeneity is markedly different from crossovers be- 
tween the three basic symmetry classes in quantum dots, 
where the functional forms are universal and given by 
random-matrix theory™ For homogeneous wires, 5G is a 
universal function of L/l so and L/Ih- 

The main assumption underlying our calculations is 
that the wire width W < L, i.e. quasi one- 
dimensionality. We also assume that the wire is well in 
the diffusive regime, I <C L,I so ,Ih <C £, where I is the 
elastic mean free path, and, for a non-homogeneous mi- 
crobridge, that the number of propagating channels at 
the Fermi level N has only one minimum along the wire 
(excluding the possibility of a "cavity"). We first dis- 
cuss our calculations for homogeneous wires; the case of 
non-homogeneous samples is discussed at the end of this 
paper. 

Starting point of our calculation is a random-matrix 
model similar to that used by Dorokhovi^ A disordered 
wire with N propagating channels at the Fermi level is 
modeled by N one-dimensional channels and periodically 
inserted scatterers that scatter within and between the 
channels. The electronic wavefunction is represented by 
a 2./V-component vector of spinors. The 2A^ components 
of the wavefunction refer to the transverse channel and to 
the left /right mover index. Linearizing the kinetic energy 
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in each of the channels, the Hamiltonian H takes the form 
of a differential operator with respect to the coordinate x 
along the wire and a 2iV-dimensional quaternion matrix 
with respect to the channel and left/right mover indices 
and spinor degree of freedom, 

d 

H = -ia (g>T 3 ®l N — + ^2v j 5{x-ja), (2) 

j 

with Co the 2x2 unit matrix for the spinor degree of free- 
dom, r 3 the Pauli matrix in left-mover/right-mover grad- 
ing, Ijy the TV x TV unit matrix in the channel grading, 
Vj a Hermitian 27V x 27V quaternion matrix representing 
the jth scatterer along the wire, and a the distance be- 
tween scatterers. A quaternion is a 2 x 2 matrix acting 
in the spinor grading with special rules for transposition 
and complex conjugation 2 ^: The "dual" X K of a quater- 
nion matrix is A R = (J2X T <72', the quaternion complex 
conjugate is defined as X* = (7V^) R . We have chosen 
units such that the Fermi velocity is one. A model sim- 
ilar to Eq. J!5J has been used in Ref. to study weak 
localization in unconventional superconducting wires. 

The ensemble-averaged conductance (G) of the wire is 
given by the Landauer formula, 



( Q ) = J9, ff=<tr(l-rtr)>, 



(3) 



where r is the TV x TV quaternion reflection matrix of 
the wire. To calculate r, we start from a wire of zero 
length and add slices of length a at the wire's ends. The 
scattering matrix of the jth scatterer is 



/ 

3 '3 



U r 
r 



2i + Vj 



(4) 



Hence, if a scatterer is added at the lead end of the wire, 
the new reflection matrix of the wire is calculated accord- 
ing to the composition rule 



rj+t'Ml-r'.rrhj 



(5) 



(A similar composition rule, involving both transmission 
and reflection matrices of the disordered wire, applies if 
a scatterer is added at the far end of the wireA) 

In left-mover/right-mover grading, the potential Vj is 
parameterized as 



V = 



VLL V LR 
VRL V RR 



(6) 



where Vlr, v R l, and v RR are TV x TV quaternion 

matrices, 



v L L{otf,r)f) = v RR (a f ,-r) f ) 



(7a) 



a 

TJn 



(u°f + rjfXf) ® (To + iaf 22 u / ® °> 



VLR(a b ,i] b ) = v RL (a b ,T] b ) 



(7b) 



l(N + l) 



1° + r) b x b ) ® cr + ia b ^ ® a^, 



In Eq. (0 , u° and x / are random Hermitian TV x TV matri- 
ces, Up fi = 1,2,3, is a random anti-Hermitian matrix, 
is a random symmetric matrix, and u%, [i = 1,2,3 
and x b are random antisymmetric matrices. All of these 
random matrices have independent and Gaussian distri- 
butions with zero mean and unit variance. (Variances 
are specified for the off-diagonal elements; diagonal ele- 
ments have double variance for symmetric matrices and 
are zero for antisymmetric matrices.) The parameters 
a b and a t describe the strength of the breaking of spin- 
rotational symmetry; The parameters rj b and r/f describe 
the strength of the breaking of time-reversal symmetry. 
Finally, I f is the elastic mean free path for forward scat- 
tering and I is the transport mean free path. 

To find the conductance of the wire we calculate the 
change of g if one scatterer is added to the wire. To this 
end, we expand the scattering matrix Sj of Eq. Q in 
powers of Vj , use the composition rule J5J , and calculate 
the Gaussian average over the potential Vj. In the limit 
a <C / of weak disorder we thus find 
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- 2Nl eL 9 = g ' 



(8) 



(9a) 



ho + 3hi. 

We abbreviated 

ho = (tr(l-rV)(l -r*r R )), 
1 3 

hi = -^<tr(l-rVK(l-r*r R K), (9b) 

and omitted terms that vanish in the diffusive regime 
I <C L,I so ,Ih <§C Nl. The subscripts and 1 refer to 
singlet and triplet contributions, respectively. 

To leading order in TV, Eq. ((SJ can be solved without 
the interference corrections ho and h\, with the result 
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2NI 

— + 0(1), 



(10) 



corresponding to the Drude law for the conductance. The 
0(1) correction in Eq. (|10f> gives the weak localization 
correction 6g, which we now compute. 

To find the weak localization correction, we need to 
calculate ho and h%. Proceeding as before, we find that 
the L-dependence of h m , m = 0, 1 is governed by the 
evolution equation 

,dh m (2NI , \ , 8TV 2 Z 2 

m ~bT = [ir + m ) m + ^~' ^ = o,i, 

(11) 

where we abbreviated 



1 3 

■o = (tr (1 - r*r)), k x = - £(tr (1 - rV„w M )). (12) 



Evolution equations for ko and k\ are obtained similarly 
and read 



2NI 



dL 

dki 



> 



(13a) 



(13b) 
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where the length scales Ih and l' H are defined in terms of 
the parameters of the random-matrix model (J7J , 



0.4 



l H 2 = 2(l- 2 ri+l-H- 1 i 1 }) 1 

I- 2 = eii-^l + rHj 1 ^), 

Cif) 2 = ^H + |^so 2 - 

Equations and l|13() have the solution 

2NI , L 
fco = - — coth—, 



h 



2NI 



l 



coth - — 
L l H 



coth 



H 



(14a) 
(14b) 
(14c) 



(15a) 
(15b) 



Expressions for k\ and h\ are obtained from Eq. i|15|) 
after the substitution Ih — > Substitution of /io and 
/ii into Eq. (JSJ then allows for the calculation of the weak- 
localization correction to the conductance, 



At zero magnetic field, Eq. IjKil) simplifies to 



1 9Z 2 

r 1 so 

53= 3 + 4^ 



21 



■ coth 



2L 



(16) 



(17) 



Equation (|17fl reproduces the limits <5G = —2e 2 /3h 
without spin-orbit scattering and SG — e 2 /ih with 
strong spin-orbit scattering. Without spin-orbit scatter- 
ing, Eq. agrees with the weak localization correction 
calculated in Ref. I22T For large magnetic fields, L^§> Ih 7 
Eq. i|lti|) simplifies to 



8g 



4 r 2 
3 so 



-1/2 



(18) 



which has the same functional form as the weak lo- 
calization obtained using diagrammatic perturbation 
theory. 14 i 15 i 16 i 23 Comparison of Eq. JTBJ) and Refs- ITITI 
Hal23l allows us to identify l so as the spin-orbit length, 
and, for a channel (with width W 3> 1) in a two- 
dimensional electron gas in a perpendicular magnetic 
field B, 



l 2 H = 3(h/WBe) 



(19) 



The case of a cylindrical wire of radius R ^> I and mag- 
netic field perpendicular to the wire is obtained by the 
substitution W 2 -> 3i? 2 /2. For I > W (or I > R) the 
crossover length Ih has a more complicated /-dependent 
expression^ 

Figure^a) shows 5g as a function of the magnetic field 
for several values of the spin-orbit coupling. In Fig. Ush) 
we show Sg as a function of l^L for several values of the 
magnetic field. 

We now turn to a description of the weak localization 
correction in a non-homogeneous microbridge. Examples 
of non-homogeneous microbridges with varying widths 
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FIG. 1: The weak localization correction 8g plotted (a) as 
a function of the magnetic field strength (characterized by 
the dimensionless ratio IJ^L) for fixed value of the spin-orbit 
scattering rate (characterized by l^ 1 L). From bottom to top, 
the curves correspond to L/l so = 0.1, 2, 4, 6, 10, 30, and 00. 
(b) as a function of length L for fixed l~^l so . From bottom to 
top, the curves correspond to lT}l so = 2, 0.3, 0.2, 0.1, and 0. 
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FIG. 2: The weak localization correction Sg as a function of 
the magnetic field strength for three different shapes of a dis- 
ordered microbridge (channels in a two-dimensional electron 
gas). The three different shapes are characterized by s(x) — 1, 
s(x) = l+4|2x/L| and s(x) = l+4(2x/L) 2 , -L/2 < x < L/2, 
cf. Eq. 121II . as shown in the inset. The three groups of curves 
correspond to strong, intermediate and weak spin-orbit scat- 
tering from top to bottom, with l so in the intermediate case 



chosen for each case to render the same correction as I , 



0. 



The magnetic field strength is measured in terms of the effec- 
tive magnetic length lij.ett, cf. Eq. 112MB . 



are shown in the inset of Fig-El If t ne w i re cross section or 
the electron density vary with the coordinate x along the 
wire, the number of propagating channels at the Fermi 
level N also varies with x. We assume that N(x) has a 
minimum for x = and that dN/dx > (dN/dx < 0) for 
all x > (x < 0). Further, x-dependence of the impu- 
rity concentration, the smoothness of the boundary, the 
shape of the cross section, etc., causes an cc-dependence 
of the length scales I, Ih, and l so . 

The reflection matrix of the wire is constructed by 
building the wire from thin slices, starting at the nar- 
rowest point x = 0. This way, the number of channels in 
the slices added to both ends of the wire can increase, but 
not decrease. For the construction of an evolution equa- 
tion for the conductance g and for the auxiliary functions 
ho, hi, fco, and fci, we distinguish between two types of 
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added slices: A thin slice that contains a scattering site 
but for which the number of channels remains constant, 
and a thin slice without scatterer in which N increases 
by unity. Addition of a slice of the former type causes 
a small change in the reflection matrix r, which is the 
same as for a quantum wire of constant thickness, see 
Eq. (JjjJ above. Addition of a slice for which N increases 
by unity does not cause a change of the conductance g 
or of the auxiliary functions ho, hi, ko, or ki, as can seen 
by inspecting the cases x > and x < separately: For 
x > 0, an increase of N does not cause additional reflec- 
tion, and hence does not affect the reflection matrix r; 
for x < 0, an increment of N changes the dimension of 
the reflection matrix r by one, 

'-(£!)■ < 2 »> 

but does not change the conductance g, or the functions 
ho, hi, ko, or k\. Combining the two types of slices, we 
conclude that the only effect of the ^-dependence of N 
and I is indirect, through the explicit appearance of N 
and I in statistics of the scattering matrix of the added 
slice, see Eq. {7J|. In the diffusive regime, N(x) and l(x) 
only appear in the combination 

s{x) = N{x)l(x)/N l , (21) 

where Nq and lo are number of propagating channels and 
mean free path at x = 0. For large N the function s(x) 
may be considered continuous, and the evolution equa- 
tions become differential equations which now include ex- 
plicit reference to the function s(x). If the wire length L 
is replaced by the effective length L, 

J s(x) 

the evolution equations for <?, ho, hi, ko, and k± keep the 
same form as for homogeneous wires, provided we make 
the substitutions N — ► No, L — + L, I — > lo, Ih — > Ih = 
l H /s(x), and l so -> l so = l so /s(x). 

The functional form of the leading-in-A contribu- 
tion to the conductance remains unchanged, G — 
(e 2 /K)(2NqIq/ L). Also, for the limiting cases of no 
spin-orbit scattering and strong spin-orbit scattering, 
the weak localization correction SG is still given by the 
universal result Eq. However, because of the in- 

dependence of the length scales Ih and l so , Sg acquires 
an explicit dependence on the shape of the disordered mi- 
crobridge or the non-homogeneity of the mean free path 
or the electron density in the crossover region between the 



symmetry classes. For a large magnetic field (l^L ^> 1), 
the weak-localization correction can be found in closed 
form, 

s 9 = 2 ( 1h ' cS ~ 31 'h,cs) i ( 23 ) 
1 f lu(x)dx , If l' H (x)dx 

Equation (|23l) simplifies to Eq. (|18l) in the case of s(x) 
constant. The same result follows if Eq. I|18() is inter- 
preted as a quantum interference correction to the one- 
dimensional resistivity and Ijj is taken ^-dependent. For 
weaker magnetic fields with l H L of order unity, a nu- 
merical solution of the evolution equations is required. 

In Fig. we show results of a numerical solution of 
5g for the examples s(x) constant, s(x) = 1 + 4\2x/L\ 
and s(x) = 1 + 4(2x/L) 2 , -L/2 < x < L/2. These 
functional forms correspond to diffusive microbridges in 
a two-dimensional electron gas of the form shown in the 
inset of Fig. [21 with uniform impurity concentration and 
mean free path I -C W. The three sets of curves in 
the figure represent strong, intermediate and weak spin- 
orbit scattering, respectively. For the intermediate case 
(middle set of curves in Fig. |2J, three different values of 
l so were chosen so that the weak-localization correction 
Sg = is equal in the three cases for zero magnetic field. 
The magnetic field is characterized by the ratio l H \ s L, 
cf. Eq. m order to remove a spurious shape depen- 
dence for the large-field asymptotes. While there is no 
dependence on the form of the function s(x) in the lim- 
iting cases of zero and large magnetic fields, we observe 
that, indeed, Sg depends on the precise form of the non- 
homogeneity for intermediate magnetic field strengths, 
although, with proper scaling, the difference between the 
results for the three cases we considered is less than 10%. 

In conclusion, we have shown that the scattering ma- 
trix approach to quasi one-dimensional weak localiza- 
tion can be used to obtain a detailed description of the 
crossover between the different universality classes. We 
have recovered some results known from diagrammatic 
perturbation theory, and have discovered one aspect of 
the problem that has not been noticed previously: The 
dependence of the functional form of the crossover on 
non-homogeneities in the conductor. 
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